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The equation of state of H2 adsorbed in the interstitial channels of a carbon nanotube bundle 
has been calculated using the diffusion Monte Carlo method. The possibility of a lattice dilation, 
induced by H2 adsorption, has been analyzed by modeling the cohesion energy of the bundle. The 
influence of factors like the interatomic potentials, the nanotube radius and the geometry of the 
channel on the bundle swelling is systematically analyzed. The most critical input is proved to be 
the C-H2 potential. Using the same model than in planar graphite, which is expected to be also 
accurate in nanotubes, the dilation is observed to be smaller than in previous estimations or even 
inexistent. H2 is highly unidimensional near the equilibrium density, the radial degree of freedom 
appearing progressively at higher densities. 

PACS numbers: 61.46.+W, 68.43.Fg, 02.70.Ss 
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I. INTRODUCTION 

Theoretical and experimental research in carbon nan- 
otubes (CN's) is nowadays a very active field in which 
new avenues are opening continuously, all related to its 
particular structure at the nanometer-scale level. Besides 
their singular electrical and mechanical properties, CN's 
are able to adsorb some atoms and molecules offering a 
real possibility for quasi-one-dimensional environments. 
These are mainly the inner part of the nanotubes and the 
interstitial channels (IC's) among them due to their nat- 
ural bundle arrangement. When CN's are formed, they 
appear as long capped cylinders which only can adsorb 
a given substance if their caps are removed by chemi- 
cal means. Instead, IC's are always present offering very 
narrow channels, distributed according to an hexagonal 
lattice, and with a sizeable adsorption capability in com- 
parison with adsorption in planar graphite. Recently, Ta- 
lapatra et aim studied experimentally the binding ener- 
gies of Xe, CH4, and Ne on close-ended nanotube bundles 
and concluded that none of them is adsorbed in the in- 
tersites. Probably, the size of these atoms and molecules 
is too large compared with the radius of the interchannel 
(ric — 3 A) to fit into them. However, the situation for 
lighter species like He or H2 seems different. There are 
several experiments which claim that both of them can be 
adsorbed in the IC's due to their small radiusiSiii H2 is 
certainly the most interesting system since a CN bundle 
could be one of the best options to reach the target en- 
ergy densities for a lightweight hydrogen-storage system 
usable for transportation^^ 

In a recent theoretical work, Calbi et al£ have studied 
the adsorption of H2, and other molecules and atoms, in 
the IC's of a CN bundle. Introducing in the formalism a 
possible dilation of the bundle, by means of a harmonic 
approximation, they concluded that in equilibrium the 



bundle can swell. The dilation is observed in all cases, 
with different intensity depending on the particular sys- 
tem. In the geometry there considered, corresponding to 
(10,10) tubes in the standard terminology, and when the 
gas adsorbed is H2 that swelling amounts approximately 
to 1 %. This is the result of a delicate balance between 
two competing effects when the lattice is expanded: a de- 
crease in the CN bundle cohesion energy and an increase 
in the H2-nanotube binding energy. Obviously, this pre- 
diction for H2 acquires special relevance since it would 
imply a significant increase of the CN bundle storage 
power. However, at present there is not any experimental 
confirmation of this possible adsorption-induced dilation. 
Probably, even if such effect is indeed present, its man- 
ifestation could be hardly observed due to the present 
experimental uncertainties. Raman spectroscopy seems 
a promising method for a better insight. In a recent work, 
this method was applied successfully by Williams et aim\ 
to discriminate the different adsorption sites of H2 in CN 
bundles. 

In the present work, we present a diffusion Monte Carlo 
(DMC) study of H 2 adsorbed in IC's of a (10,10) CN 
bundle. In the past, and using the same methodology, 
we characterized the ground-state properties of 4 He and 
H2 inside a single CN2*S and of 4 He adsorbed in an IC>i£ 
In this paper, our aim is to determine the equation of 
state of H2 using microscopic theory, with special atten- 
tion to the possibility of an adsorption-induced bundle 
dilation. To this end, the influence on swelling of in- 
puts like the nanotube radius, the molecule-molecule and 
CN-molecule interactions, the geometry of the intersite, 
and the transversal degree of freedom are thoroughly dis- 
cussed. Our results agree somehow with the ones of Ref. 
7 but the magnitude of a possible swelling is seen to be 
comparatively smaller or even inexistent. 

The rest of the paper is organized as follows. In Sec- 
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tion II, the influence of the several parameters of the 
model on a possible dilation of the bundle is studied us- 
ing a one-dimensional approximation. The accuracy of 
this model is then tested in Section III by means of a 
three-dimensional DMC calculation. Finally, Section IV 
comprises a brief summary and the main conclusions. 



II. ONE-DIMENSIONAL APPROXIMATION: 
TESTING THE MODEL 

As in Ref. 0, we study H2 adsorbed in the IC of a 
(10,10) CN bundle. The bundle is disposed in such a way 
that, in the plane perpendicular to the CN long axis, a 
triangular lattice is formed. The lattice constant, i.e., 
the distance between adjacent centers, is Do = 17 A. Do 
corresponds to the equilibrium position and there is an 
overall agreement on its value. Therefore, we consider it 
as a fixed parameter in all the calculations. 

Having established the value of Do, the radius of a 
single nanotube is directly related to the C-C minimum 
separation between two neighbors. In a theoretical calcu- 
lation, Tersoff and Ruoff 11 concluded that this distance is 
3.4 A, in agreement with the most accepted value for the 
Lennard- Jones uc-c- In fact, this value coincides with 
the experimentally accepted distance between graphite 
sheets. Considering this C-C distance, the radius of a 
(10,10) tube is 6.8 A. It is worth mentioning that this 
radius is the same one obtains constructing the same 
nanotube by rolling up a graphitic surface with a C-C 
distance of 1.42 A. However, there is not a general agree- 
ment on those values. Recently, Charlier et alMi carried 
out a density functional calculation of a CN bundle and 
obtained an equilibrium geometry corresponding to a C- 
C intertube distance of ~ 3.2 A. The CN radius is then 
6.9 A, a value which was used in Ref. 0. At present 
and to our knowledge, there is not experimental data on 
the C-C distance which could help to choose between dif- 
ferent theoretical models. Therefore, we consider in the 
present calculation two possible CN radius, 6.8 and 6.9 

A. 

The interstitial channel between three adjacent nan- 
otubes in a (10,10) bundle can lodge a hard sphere with 
radius 2.9 or 3 A, corresponding to CN's with radius 6.8 
and 6.9 A, respectively. Since the a parameter for the C- 
H2 Lennard- Jones potential is around 3 A, one can rea- 
sonably assume a one-dimensional (ID) approximation 
for hydrogen adsorption in IC's. In this approximation, 
the H2 energy per particle (e = E/N) can be written as 7 - 

e(A, D) = e (D) + e 1D (A) + h(X, D) + j^k(D~ D ) 2 , 

(1) 

with A the linear density of the H2 molecules and D the 
lattice parameter. In Eq. (JIJ, £o(D) is the binding energy 
of a single molecule in the IC, eiD(A) is the ID H 2 energy 
per particle, /i(A, D) corresponds to the interaction en- 
ergy with H2 in neighboring IC's, and the last term takes 
into account the cohesion energy of the bundle when it is 



dilated from the equilibrium lattice distance -Do- In the 
following, we analyze the several terms entering Eq. Q 
to disentangle which are the relevant inputs influencing 
a possible dilation of the bundle. 

Both e (D) and eiD(A) are calculated using the diffu- 
sion Monte Carlo (DMC) methodiiii The DMC method 
solves in a stochastic way the V-body Schrodinger equa- 
tion in imaginary time for the wave function f(R,t) = 
1>(R)*(R,t), 



df(R,t) 
dt 



= -DV 2 f(R,t)+DV(Ff(R,t)) (2) 

+(MR)-E)f(R,t) , 



with ip{R) a trial wave function introduced for im- 
portance sampling. In Eq. (J2J, D = fi 2 /2m, 
El(R) = ijj(R)^ 1 H^p(R) is the local energy, and 
F(R) = 2ip(R)- 1 Vip(R) is the so-called drift force 
which guides the diffusion movement to regions where 
i/j(R) is large. 

The trial wave function ip used in the calculation of ID 
H2 is a Jastrow one, 



ip 1B (R) = Hf(n J ) 



(3) 



i<3 



with f(rij) = exp (-0.5 (&/%) J a McMillan two-body 
correlation factor. The parameter b in /(r) has been de- 
termined by means of a variational Monte Carlo (VMC) 
optimization. Near the equilibrium density, the optimal 
value is b = 3.996 A and it increases gradually with the 
density A (at A = 0.358 A -1 , b = 4.026 A). 

The trial wave function used in the calculation of £o(D) 
includes two-body correlations with the three nanotubes 
surrounding the IC. They are also of McMillan type, 



^ ic (R)= n^xp 



(4) 



r n being the distance of the hydrogen molecule to the cen- 
ter of any of the three tubes. The parameter a, optimized 
using VMC, varies from a = 22.5 A for zero dilation to 
a = 21.5 A for D - D Q = 0.2 A. On the other hand, the 
optimal values for the parameter a present a negligible 
dependence with the particular C-H2 potential chosen in 
the calculation. The IC is so narrow that one can ap- 
proximate the sum of individual CN-molecule potentials 
by a new one which is only a function of the radial dis- 
tance r to the center of the IC. This simplified model is 
obtained by an azimuthal average of the three individual 
potentials. 7 In this case, the trial wave function is simpler 
than the previous model @. We have chosen a gaussian 



tp c (R) = exp(-cr2) , 



(5) 



the parameter c varying from 7.6 A -2 to 11.5 A~ 2 , de- 
pending on the radius of the tube and on the dilation. 
The greater the radius, the smaller the value of c. 
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Para- hydrogen is spherical to a large extent. As usual 
in most of molecular hydrogen calculations, the inter- 
molecular interaction is considered purely radial and 
described by the Silvera and Goldman (hereafter SG) 
model. 15 In addition, and to make comparisons with pre- 
vious work (Ref. 0), we have also made some calcula- 
tions with the potential proposed by Kostov, Cole, Lewis, 
Diep, and Johnson>i£ This H2-H2 potential (hereafter 
KCLDJ) incorporates, in a rather crude way, three-body 
corrections to the pair potential coming from the triplets 
H 2 -C-H 2 . 

Much more critical for estimating a possible bundle di- 
lation is the model for the CN-H2 interaction. As in our 
previous work on the equation of state of H2 adsorbed in- 
side a nanotubejS we use the CN-H2 potential proposed 
by Stan and Cole>i£ This potential results from an av- 
erage over all the C-H2 interactions between the C of 
an infinite CN and a single molecule located at a radial 
distance r. The interaction is thereby independent of 
z (corrugation effects are neglected) and can be used for 
adsorption both inside and outside the nanotube. The re- 
sulting CN-H2 potential depends explicitly on the a and 
e parameters of the Lennard-Jones potential between C 
and H2. At present, there is not an overall agreement 
about which are the best set of parameters (cr,e) describ- 
ing this interaction. In order to study the influence of 
this choice on the calculations, two different sets have 
been studied. The first one, used for example in Refs. 
andll8l is derived with the Lorentz-Berthelot combining 
rules 



£c-h 2 — V £h 2-h 2 e C-C ■ 

Using ctc-c = 3.4 A, e c _ c = 28 K, cth 2 -h 2 = 3.05 
A, and £h 2 -h 2 = 37.0 K the first set is ec-H 2 = 32.2 
K and ctc-h 2 = 3.23 A(hereafter LB). The second op- 
tion is to consider the optimal parameters which de- 
scribes the interaction of H 2 with planar graphite. This 
set from Wang, Senbetu, and Woo^ is probably more 
realistic than the first one and it has also been used 
in the pasti 3 i 9 i 20 The values are ec-H 2 = 42.8 K and 
cr c _H 2 = 2.97 A(hereafter WSW). Notice the sizeable dif- 
ference between both sets which will generate significant 
differences in the final results. 

Dilation and compression of the lattice constant in the 
CN bundle have an energetic cost. Following Refs. |^] 
and 7, that contribution is assumed to be harmonic in 
the displacement around the equilibrium position Dq (see 
Eq. ^| . This approximation is expected to hold for small 
D — Dq values but the uncertainty in the real value of 
k is rather large. From compression modes measured 
in graphite by Nicklow et al.^ Mizel et a/iSl estimated 
a value k — 1740 K A~ 3 . The constant k can also be 
obtained from the bulk modulus B. Considering that 
any deformation of the tube takes place only in a di- 
rection perpendicular to its long axis, 
From the experimental measure of B by Tang et alm& 



Radius (A) IC V C -H 2 eg (K) 

6.9 T WSW -1020.69(8) 

6.9 T LB -278.02(8) 

6.8 T WSW -1096.02(2) 

6.8 T LB -668.39(9) 

6.9 R WSW -965.93(7) 
6.9 R LB -207.59(5) 
6.8 R WSW -1044.49(2) 
6.8 R LB -576.56(2) 

TABLE I: Binding energy of a single H2 molecule adsorbed in 
the interstitial channels of a (10,10) CN bundle for different 
CN radius, IC geometries, and C-H2 potentials. Figures in 
parenthesis are statistical errors. 

(B = 41.66 GPa) one obtains k = 3015 K A" 3 , a fig- 
ure much larger than the one from Ref. 0. Part of the 
difference can somehow emerge from the fact that the ex- 
perimental data for B was obtained for CN's with wider 
radius (7.04 A). In the results presented below, we have 
used k = 1740 K A -3 everywhere but the influence on 
the final results of the particular k value is analyzed in 
selected cases. 

Finally, the term h(X, D) in Eq. Q sums the interac- 
tion energy between the H2 molecules of a given IC and 
the ones of its neighboring channels. This contribution 
may be readily obtained assuming a mean field approxi- 
mation, 

\ r°° / N 

h(X,D) = - J dx Uh 2 -h 2 (Vx 2 + £ 2 ) , (7) 

which considers the neighboring channels as uniform ar- 
rays without correlation effects that could modify the H2 
kinetic energy. The calculation in Eq. Q is extended to 
the nearest neighbors, the next-nearest neighbors, and 
then on up to the desired accuracy. Obviously, h(\,D) 
changes when the bundle swells. In Ref. |2J, we have 
verified the high accuracy of Eq. Q by comparing its 
estimation with an exact DMC calculation. 

The influence of the CN radius and the C-H2 poten- 
tial in the equation of state of adsorbed H2 is clearly 
observed in the value of eo- Table I contains results for 
the H2 binding energy considering two CN radius, 6.8 
and 6.9 A, the two C-H2 potentials discussed above, LB 
and WSW, and two geometries for the IC, the real one 
with the triangular composition of the tubes surrounding 
the IC (T) and the radial (R) model from an azhimu- 
tal average. The narrowness of the IC makes the differ- 
ences between the T and R results not to be larger than 
15 %. The dominant effect is unquestionably the C-H2 
potential. The significant differences between the pairs 
(<7,e) for the two models generate a big discrepancy in 
the value of eo. The largest binding energies correspond 
to the most accurate WSW potential due to its larger 
parameter e. It is also noticeable, and will be discussed 
below in connection with dilation, the different behavior 
of e with the CN radius: it strongly depends on the ra- 
dius for LB whereas remains nearly unchanged for WSW. 
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FIG. 1: Energy per hydrogen molecule for CN radius 6.9 
A, T geometry, SG H2-H2 potential, and LB C-H2 potential. 
iFrom top to bottom in A = 0.2 A -1 : D - D = 0.18, 0.17, 
0.16 and 0.15 A. 
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FIG. 2: Same than in Fig. 1 but changing the LB C-H2 
potential by the WSW one. From top to bottom at A = 0.2 
A" 1 , lattice dilations of 0.06, 0.05 and 0.04 A. 



This last feature is a consequence of the larger a value 
of the LB potential which effectively reduces the space 
available inside the IC to accommodate the H2 molecule. 

In order to determine if a dilation of the bundle is ener- 
getically preferred, a series of calculations using Eq. ||TJ 
has been carried out. As commented before, our aim is to 
quantitatively determine the influence of the potentials, 
the CN radius, the elastic constant k, and the geometry 
model for the IC on a possible dilation. The direct output 
of Eq. Q is the energy per particle as a function of the 
linear density A and D. For a given dilation D — Do > 0, 
with fixed D = 17 A, e(A) can be obtained and from it 
the equilibrium point (Ao, erj) corresponding to zero pres- 
sure. Illustrative outputs of this procedure are shown in 
Figs. 1 and 2. Fig. 1 contains the H2 equation of state 
in the IC for values of the lattice dilation D — Dq = 
0.18, 0.17, 0.16, and 0.15 A. All the curves have been 
calculated assuming a CN radius of 6.9 A, the SG H 2 -H 2 
potential, k = 1740 K A~ 3 , the LB C-H 2 potential, and 
a T channel. The lowest energy per particle at the mini- 
mum is obtained for a dilation 0.16 A, nearly coincident 
with the result reported in Ref. T (0.166 A) obtained 
with the H2-H2 KCLDJ potential and an R channel. If 
the C-H 2 potential is the WSW model, the H 2 equation 
of state changes dramatically and the dilation becomes 
much smaller. This is shown in Fig. 2 which differs from 
Fig. 1 in the C-H 2 potential, the other inputs being un- 
changed. The curves now correspond to dilations D — Dq 



— 0.06, 0.05, and 0.04 A. The minimum corresponds to a 
swelling of 0.05 A. The equilibrium density is A = 0.232 
A~\ slightly inferior to the one obtained with the LB 
potential, A = 0.268 A -1 . The magnitude of the energy 
per particle is very different in the two cases: going from 
LB to WSW the energy increases in a factor two. The 
major part of this increment comes from the differences 
in the binding energies of a single molecule eo (Table I). 

The full set of results combining CN radius, C-H 2 and 
H 2 -H 2 interactions, and IC geometry is reported in Ta- 
ble II. Ao corresponds to the H 2 equilibrium density at 
the energetically preferred swelling, reported also in the 
Table. Due to the large uncertainty in the value of k, 
we only report the influence of its change (3015 K A~ 3 
instead of 1740 K A™ 3 ) in the most favorable cases for 
the bundle swelling. Inspection of the data shows that 
neither the H 2 -H 2 potential nor the detailed form of the 
IC, R or T, are relevant in the swelling process. Main- 
taining fixed the radius to a value 6.9 A and the LB 
C-H 2 interaction, the dilation comes out 0.16-0.17 A, in 
overall agreement with Calbi et al.J However, if the LB 
interaction is substituted by the WSW one, the swelling 
observed is reduced by a factor of three (from 0.17 A to 
0.05 A), and the equilibrium density Ao decreases by 15 
%. 

The second part of Table II contains a similar analysis 
for a CN radius of 6.8 A. A general trend arising from 
the comparison between the first half of the Table (6.9 A) 
and the second one (6.8 A) is the increase of the dilation 
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Radius (A) IC V C -h 2 Vh 2 -h 2 (D 


- Do) 


(A) A (A 


eo (K) 


6.9 


T 


WSW KCLDJ 


0.04 


0.2087(3) 


-1032.98(1) 


6.9 


T 


LB KCLDJ 


0.16 


0.2667(1) 


-481.42(2) 


6.9 


T 


WSW SG 


0.05 


0.2322(2) 


-1041.41(1) 


6.9 


T 


LB SG 


0.17 


0.2685(1) 


-491.64(1) 


6.9* 


T 


LB SG 


0.11 


0.276(12) 


-300.35(1) 


6.9 


R 


WSW KCLDJ 


0.05 


0.2252(1) 


-980.07(1) 


6.9 


R 


LB KCLDJ 


0.16 


0.2662(1) 


-422.29(2) 


6.9 


R 


WSW SG 


0.05 


0.2323(3) 


-988.59(1) 


6.9 


R 


LB SG 


0.16 


0.2685(1) 


-432.17(1) 


6.8 


T 


WSW KCLDJ 


0.01 


0.1311(6) 


-1096.35(1) 


6.8 


T 


LB KCLDJ 


0.10 


0.2472(1) 


-731.55(2) 


6.8+ 


T 


WSW SG 


0.01 


0.2194(5) 


-1102.31(2) 


6.8 


T 


LB SG 


0.01 


0.2514(2) 


-740.64(2) 


6.8 


R 


WSW KCLDJ 


0.01 


0.1315(2) 


-1045.85(1) 


6.8 


R 


LB KCLDJ 


0.11 


0.2501(5) 


-679.27(1) 


6.8 


R 


WSW SG 


0.02 


0.223(1) 


-1047.39(2) 


6.8 


R 


LB SG 


0.10 


0.2502(9) 


-679.53(3) 



TABLE II: Dilation ((D — Do)), equilibrium densities (Ao) and energies eo for different radius, H2-H2 and C-H2 interactions, 
and geometries (T and R) . Figures in parenthesis are statistical errors. 
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FIG. 3: Full DMC calculation of molecular hydrogen ad- 
sorbed in the interchannels of a CN bundle. No dilation, full 
squares; D — Do — 0.05 A, empty squares; D — Do = 0.10 A, 
full circles; D — Do = 0.15 A, empty circles. 



with the CN radius when the rest of the inputs are kept 
fixed. This seems quite obvious since, with a smaller 
nanotube and the same CN-CN distance, the room for 
adsorption increases. Thus, the swelling decreases from 
0.17 to 0.11 A with the LB C-H2 interaction whereas it 
turns absolutely negligible in case of using the probably 
most accurate WSW potential. On the other hand, Aq 



is observed to systematically decrease when the radius 
is moved from 6.9 to 6.8 A. Finally, the influence of the 
value of k appears in the Table marked with an aster- 
isk for a particular case in which the dilation is one of 
the largest values. As it could be expected, (D — D ) 
decreases if k increases: using k = 3015 K A~ 3 instead 
of 1740 K A~ 3 the dilation is reduced by 30 %. Similar 
reductions would be obtained in the other cases. 



III. FULL DMC CALCULATION 

In the preceding Section it has been assumed that ra- 
dial and longitudinal degrees of freedom of H2 inside IC's 
are not coupled. The correlations between H2 molecules 
were purely longitudinal and the interaction with the sur- 
rounding walls was solved only for the one-body problem. 
In this Section, we check the validity of that approxima- 
tion by making an exact three-dimensional DMC calcula- 
tion of H2 in the IC. The CN bundle cohesion term and 
the mean field contribution h(X, D) are summed up to 
the DMC energy as in Eq. (|TJ. 

The trial wave function for importance sampling is 
written as 



i>(R) = i> w (R)^ c (R 



(8) 



where -0 1D (i?) and ip (R) are the same than in the ID 
calculation, Eqs. © and 10)), respectively. The optimal 
values for the parameters entering Eq. JSJl are the same 
than the ones in the preceding Section. The DMC cal- 
culation has been carried out for a selected case which 
we consider contains the most reasonable model. The 
CN radius is 6.8 A and the IC is of T type. The H 2 -H 2 
interaction is SG and the parameter set for C-H2 is the 
WSW one. 

In Fig. 3, the total energy per particle of H2 in the 
IC is shown as a function of A and for different bundle 



6 



10 n 1 1 1 1 1 1 r 

8 - 

6 - 

^ 4 - 

>> o / ■ 

bfj 2 - 



-2 - 

■ 

/ ■ 

-4 - 

-6 ^ 1 " ' * ' 1 1 1 1 1 

0.2 0.22 0.24 0.26 0.28 

density (A -1 ) 

FIG. 4: Comparison between a DMC calculation (full 
squares) and a ID approximation (solid line) in a case without 
dilation. See comments in the text. 

lattice values, D - D = 0, 0.05, 0.10, and 0.15 A. The 
minimum of the energy is achieved with D = Dq, what 
implies an absence of dilation in the bundle of nanotubes. 
This is in agreement with the corresponding case in Ta- 
ble II, marked by a dagger. In the full 3D calculation 
the equilibrium density is A = 0.2184(12) A^ 1 , to be 
compared with 0.2194(5) A -1 from Table II. The second 
value is within the error bar of the first one. The en- 
ergy per hydrogen molecule at the equilibrium density is 
also very similar in both cases, eo = —1101.43(1) K and 
eJ D = —1102.31(2) K. This clear overlap means that the 
ID model is a very good approximation to the physics 
of H2 inside the intcrchannels of a bundle of nanotubes 
around the equilibrium density. 

Finally, in Fig. 4 the comparison between both calcula- 
tions is extended to higher densities. In the energy scale, 
the binding energy eo = —1096.02 K (Table I) is now 
subtracted. Equation (JIJ provides un upper bound to 
the exact result which is very close to the exact energies 
near Ao- However, its quality worsens when the density 
increases due to the emergence of the radial degree of 
freedom which is frozen in the ID model. 



IV. SUMMARY AND CONCLUSIONS 

The equation of state of H2 adsorbed in the IC's of a 
CN bundle has been calculated using microscopic theory. 



By means of the DMC method, which provides exact re- 
sults for given model potentials, the possible dilation of 
the bundle induced by adsorption has been carefully an- 
alyzed. The most complete analysis has been done in 
the ID modeli In spite of its simplicity, this approach 
has proven to be very accurate when compared with a 
full 3D calculation, specially near the equilibrium density. 
Playing with the different alternatives for the parame- 
ters entering into the calculation the influence of each 
one has been established. Summarizing, the results show 
that the critical one is the particular C-H2 interaction. 
The deeper and probably more realistic WSW potential 
reduces the swelling predicted by the LB model in a sig- 
nificant way. The same can be concluded about the C-C 
interaction, represented here by the k parameter, but to 
a lesser extent. On the other hand, sizeable changes in 
the H2-H2 interaction do not modify the physics of the 
dilation process. Our results do not exclude completely 
the possibility of some lattice dilation but if it exists its 
value seems clearly smaller than the figures quoted in 
Ref. 

A key point in the discussion is the issue of what C-H2 
potential is more realistic. To this end, one should com- 
pare theoretical data with experimental results. How- 
ever, the latter are really scarce. Recently, Vilches and 
collaborator9^i2& reported data on H2 isotherms on bun- 
dles of closed-capped carbon nanotubes for several tem- 
peratures. At very low coverages, they deduced the exis- 
tence of H 2 adsorption on the ridges created between two 
nanotubes in the outer part of the bundle and/or in the 
IC's. For these sites, they estimated a binding energy ~ 
700 K, 1.5 times larger than the one for molecular hy- 
drogen on graphite. Comparison with Table I rules out 
the case with R= 6.9 A and the LB potential: its eo is 
too small. On the other hand, the value for R = 6.8 
A and the LB potential seems to fit perfectly well the 
experimental result. However, there is a problem: the 
binding energy 700 K is supposed to be an average of 
the energies of hydrogen adsorbed on the ridges and in 
the IC's. Then, the real binding energy in IC's should be 
greater than the experimental finding, fitting well with 
the results obtained with the WSW potential. 
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